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The existence of a fundamental ultraviolet scale, such as the Planck scale, may lead to modifications 
of the dispersion relations for particles at high energies, in some scenarios of quantum gravity. We 
apply effective field theory to this problem and identify dimension 5 operators that do not mix with 
dimensions 3 and 4 and lead to cubic modifications of dispersion relations for scalars, fermions, and 
vector particles. Further we show that, for electrons, photons and light quarks, clock comparison 
experiments bound these operators at 10 _5 /Mpi. 



I. INTRODUCTION 

The Planck mass, Mpj, a dimensional parameter deter- 
mining the strength of gravitational interaction, remains 
a source of conceptual problems for quantum field theo- 
ries. When the momentum transfer in a process of two 
particle collision is comparable to the Planck mass, the 
graviton exchange becomes strong, signifying a break- 
down of the perturbative field theory description. Even 
without having a fully consistent fundamental theory at 
hand, one can hypothesize several broad categories of 
the low-energy effects induced by Mp\. The first group 
of such theories has only minor modifications due to the 
existence of new physics at Mpj. By minor modifications 
we understand that all "sacred" symmetries (Lorentz 
symmetry, CPT, spin-statistics, etc.) of the field the- 
ory remain unbroken at low energies. String theory in 
the critical dimension and in the simplest background 
reduces to field theory and gravity at the scales lower 
than M s , and provides a perfect example of this category. 
In this case the chances to probe 1/Mpi suppressed ef- 
fects are very remote, as one has to have a unrealistically 
large momentum transfers. Consequently, the propaga- 
tion of a free particle with large energy/momentum is 
immune to the effects of new physics, as all corrections 
to the dispersion relation could be cast in the form of 
(p 2 ){p 2n /M^) = (m 2n+2 /M|f) where m is the mass of 
the particle and p is the four-momentum. 

In the second class of scenarios, that include loop quan- 
tum gravity, the Mpi effects have much more "vivid" 
properties. In this approach, the discrete nature of space 
at short distances is expected to induce violations of 
Lorentz invariance and CPT. Such violations are also of- 
ten discussed in a broader context using field theoretical 
language [1] . Here one might assume a perfectly Lorentz 
symmetric, CPT conserving action at 1/-Mpj order, and 
account for the existence of the Lorentz breaking terms 
via a set of higher dimension operators. This would lead 
to modifications of the dispersion relation for a free par- 
ticle as terms of the form E n+2 / 'Mpj can appear. Such 
effects can be searched for with astrophysical tools [2-5] 
or with high precision low-energy experiments [6,7]. 



In this Letter, we consider cubic modifications of dis- 
persion relations, which would appear at the leading or- 
der 1/Afpi and which have received considerable atten- 
tion in the literature recently [2-5]. In the language of 
the effective field theory, such modifications can be de- 
scribed by the dimension 5 operators. Although dimen- 
sion 3,4 operators were extensively studied [1], dimen- 
sion 5 operators remain unclassified and poorly explored. 
Lorentz breaking can be achieved by the introduction 
of an "external" background with open Lorentz indices. 
While more complicated choices are conceivable, we con- 
sider the case where a preferred frame is defined by a 
fixed time- like four- vector n a . Within this framework, 
we study cubic modifications of the dispersion relation 
for scalars, vectors and fermions. Specializing the op- 
erators to Standard Model particles, we show that the 
cubic modification is not possible for the Higgs parti- 
cles, must have an opposite sign for opposite chiralities 
of photons, and is independent for different chiralities of 
fermions. Note that the commonly studied dispersion 
relations [2-4] do not conform with those derived here. 

Further we consider certain experimental bounds on 
these operators. It is important to realize that the intro- 
duction of any preferred frame leads to spatial anisotropy. 
The four- vector n a that characterizes the preferred frame 
may have a dominant component n° in some cosmic 
frame, but the motion of our galaxy, Solar system and 
Earth will create spatial components n, ~ 1CP 3 for a 
terrestrial observer [8]. Hence clock comparison experi- 
ments [9] or experiments with spin-polarized matter [10] 
searching for the spatial anisotropy can impose stringent 
bounds on violations of Lorentz symmetry in this con- 
text. Recently direction-sensitive clock comparison re- 
sults were used to constrain the dispersion relation for 
nucleons [7] . Here with a direct application of the spin 
precession bounds, we derive O(l) limits for leptons, 
O(10 -5 ) for photons and O(10~ 8 ) for quarks in units 
of 1/Mpi. The results in the lepton sector can be im- 
proved to the level of 10 _5 Af^ with a renormalization 
group analysis. Comparing to recent calculations [11-13], 
we find that these bounds already present an interesting 
challenge for loop quantum gravity. 
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II. CUBIC MODIFICATIONS OF DISPERSION 
RELATIONS BY DIMENSION 5 OPERATORS 

In the framework of low energy effective field theory, 
the modified dispersion relations should be derived from 
some appropriate modification of the kinetic terms in the 
Lagrangian. In this paper we shall assume that short- 
distance physics does not generate dimension three and 
four terms directly. Even though this assumption has se- 
rious theoretical problems (e.g., one wants a symmetry to 
forbid such operators), it is quite safe from the practical 
side, because even if these terms exist they are already 
severely constrained by experiment (see, e.g., [1]). 

At the next level are dimension five operators which 
would lead to 0(E 3 ) modifications of the dispersion re- 
lations. We adopt the simplest approach where Lorentz 
symmetry is broken a background four vector n a (with 
n ■ n = 1). We construct operators satisfying six generic 
criteria: 

1. Quadratic in the same field 

2. One more derivative than the usual kinetic term 

3. Gauge invariant 

4. Lorentz invariant, except for the appearance of n° 

5. Not reducible to lower dimension operators by 
the equations of motion 

6. Not reducible to a total derivative 

Conditions 2 and 5 ensure that these operators lead to 
0(E 3 ) modifications of the dispersion relations, rather 
than 0(E 2 m) or 0(Em 2 ), where m is the mass of the 
particle. Our working assumption will also be that these 
operators are naturally suppressed by a factor of 1/Mpi, 
and that m,E <C Mp\. This scaling ensures that all 
operators of dimension five can be regarded as small per- 
turbations. We shall consider cases of a scalar, fermion 
and vector particle. 

Scalar: We work with a complex scalar field with the 
standard dimension four Lagrangian: £„ = |9$| 2 — 
m 2 |$| 2 . Hence the leading order equation of motion 
yields (□ + m 2 )<f> = 0, or in momentum space with 
$ ~ cxp(-ifc • x), (-fc 2 + m 2 )$(fc) = 0. 

It is relatively easy to see that there is only one possible 
term which appears at the dimension five level 

C s = l ^(n-d) 3 <S> (1) 
Mpi 

which is odd under CPT and charge conjugation. The 
corresponding equation of motion becomes 

(□ + m 2 )$ = z-^(n-d) 3 $. (2) 
Mpi 

Another possible operator, $n • 9n<f>, is reducible to 
m 2 §n-d<& and does not effect the dispersion relation in a 
significant way. In momentum space (with n ■ d ~ —iE), 
and in the Lorentz frame where n a = (1, 0, 0, 0), we have 
the dispersion relation 

E 2 c\p\ 2 +m 2 + ^-\p\ 3 (3) 



where we have used E ~ \p\ for high energies. 

Note that for a real scalar (like the Standard Model 
Higgs), the operator C s is the total derivative and does 
not produce a cubic modification of the dispersion rela- 
tion. For the case of complex scalar, the introduction of 
C 3 respects the phase invariance of C a . The dimension 
five operator breaking this symmetry, e.g., <!>(n • <9) 3< 3> 
vanishes again as a total derivative and is forbidden if 
the phase symmetry is gauged. 

Vector: We work with a U(l) gauge field for which 
the leading kinetic term is C = —F 2 /A. (Extending 
the following to a nonabelian vector is straightforward.) 
Hence the leading order equations of motion are just 
the Maxwell equations: d a F ab = 0. After gauge fixing 
d ■ A = 0, this yields OA a = 0, or in momentum space 
k 2 A a (k) = 0. Again, we wish to modify the dispersion 
relation at 0(E 3 ) and so the new terms should satisfy 
the constraints listed above. Keeping in mind the lead- 
ing order Maxwell equations in vacuum and the Bianchi 
identities <9[ a i<fc c ] = 0, one finds that there is only one 
term with the desired form and which produces a non- 
trivial modification of the dispersion relations 

c ~< = -^r n<1F ^ n ■ d ( n bF hd ), (4) 

Mpi 

where F ab = \e ahcd F cd . Note that this operator is odd 
under CPT and even under charge conjugation. The 
equation of motion becomes 

aA a = --^e abcd n»(n-d) 2 F cd (5) 
Mpi 

where we used Bianchi identity and d- A = 0. To identify 
the effect of the new term on the dispersion relation, we 
go to momentum space and select photons moving along 
the z axis with k a — (E,0,0,p). Then for transverse 
polarizations along the x and y axes, 

( E " - P 2 ± J^-P 3 ) (e*±ie w ) = (6) 

where we have used E ~ p to leading order and chosen 
the "rest frame" for n a as before. Hence the sign of the 
cubic term is determined by the chirality (or circular po- 
larization) of the photons. As in the case of optical activ- 
ity, this leads to the rotation of the plane of polarization 
for linearly polarized photons, which is used to bound £ 
in [5]. Given that is unique, the common approach 
[2-4] of postulating a cubic dispersion relation which is 
chirality independent is incompatible with effective field 
theory. 

Spinor: We work with a Dirac spinor for which the lead- 
ing kinetic term is: C a = ^(i^ — m)^. Again, we wish to 
modify the dispersion relation at 0(E 3 ) and so consider 
new terms satisfying the constraints listed above. Keep- 
ing in mind the leading order equations, one finds that 
there are only two terms with the desired form 
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C f = —^( m 4 + m 4j 5 )(n-d) 2 ^ (7) 

Both operators break CPT, with r/i being charge con- 
jugation odd and 772 charge conjugation even. Now the 
equation of motion takes the form 

(i$-m)V = --^-{riii + milh) (n- d) 2 ^ (8) 
Mpi 

To produce an expression that is readily identified with 
a dispersion relation, one operates on both sides with 
(i@ + m) to produce 
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(vi + mih) (n ■ a) 3 * (9) 



where we have again dropped terms of order m/Mp\. 
Hence the modified dispersion relation becomes 

E 2 \p\ 2 -m 2 - ^-(m + 77275)) * = (10) 

with E ~ |p| for high energies. At high energies (i.e., 
E 2 3> m 2 ), we can choose spinors as eigenspinors of 
the chirality operator and redefine coupling constants as 
Vr.l — Vi i V2- I n the extension of the Standard Model, 
the chiral choice for 77 couplings is prescribed by gauge 
invariance. Previous studies [3,4] consider only chiral- 
ity independent dispersion relations for fermions and so 
implicitly fix 772 = 0. 

Above we have identified interesting operators which 
modify the dispersion relations at cubic order. One could 
also consider frame-dependent modifications of the inter- 
action terms between, e.g., photons and electrons. How- 
ever, these would not effect the threshold tests which are 
usually studied, e.g., [3]. On the other hand, they may 
play a role in indirect tests as those considered below. 

As a further technical aside, one can see a very seri- 
ous potential problem with divergences at the loop level. 
For example, it is easy to show that when inserted in a 
self-energy loop, the 772 operator generates a dimension 
3 operator A^yMpj 1 V> with Aj/y, the ultra-violet 
cutoff on the momentum integral. The latter cannot be 
lower than the electroweak scale or SUSY breaking scale. 
However, even assuming such a low cutoff, the bounds 
of [1] on these dimension three operators restrict all the 
couplings, T]i,^,K to be smaller than 10~ 10 . However, 
we may evade this bound as follows: The external ten- 
sor n a n b n c appearing in all of the operators above may 
be separated into a vector piece and a traceless sym- 
metric tensor C' abc = n a n b n c - i(n a g 6c +cyclic). The 
quadratic divergences above are associated with the vec- 
tor part only, while the operators proportional to C abc 
can generate only logarithmic divergences which should 
be interpreted as the renormalization group evolution of 
these operators. Thus, to be consistent with the assump- 
tion about the absence of large dimension 3 or 4 oper- 
ators, we make the substitution n a n b n c — ► C abc in the 
following. Note that this change does not affect the dis- 
persion relations in the regime E » m. 



III. LOW-ENERGY LIMITS ON DIMENSION 
FIVE OPERATORS 

Modifications of dispersion relations for stable parti- 
cles such as electrons, light quarks, and photons could be 
searched for using the astrophysical probes [2-5] . Here we 
show that terrestrial limits imposed by the clock compar- 
ison experiments are equally and sometimes more sensi- 
tive to dimension 5 operators. Indirect limits exploit the 
idea that the choice of n a = (1, 0, 0, 0) defines a preferred 
frame that does not coincide with the laboratory frame 
on the Earth [8]. Assuming that the rest frame for n a is 
related somehow to either the cosmic or galactic frame, 
the typical size of the spatial component of n is 10 -3 . 

Limits on the operators involving electrons and elec- 
tron neutrinos are especially easy to derive. We use 
the fact that best tests of directional sensitivity in the 
precession of electrons limit the size of interaction be- 
tween the external direction and the electron spin at the 
level of 10 -28 GeV [10]. This immediately translates into 
the following limit on the coefficients i]l and i]r that 
parametrize the effective interaction of the form (7) for 
left-handed leptons and right-handed electrons: 



Wl ~ Vr 



10- 28 GeVMp! „ 
< — — ■ ~ 4, 



(11) 



where we use Mpi = 10 19 GeV. The square of the elec- 
tron mass originates from (n°d t ) 2 acting on the electron 
wave function. The orthogonal combination, ri£ + rf R , is 
almost not constrained, as it does not contribute into the 
electron spin Hamiltonian. Similarly, we do not constrain 
operators involving second and third generation leptons. 
Note that our limit is comparable to the existing con- 
straints [3] . However the latter analysis assumed 772 = 
and so our result is a complementary constraint. 

The absence of a preferred direction is checked with 
even greater precision using nuclear spin, which trans- 
lates into more stringent limits on new operators for the 
light quarks. The photon operator (4) will also contribute 
because of the electromagnetic interactions inside the nu- 
cleon. To use the best experimental limits of 10~ 31 GeV 
on the coupling of rij to neutron spins [9] , we must relate 
the photon and quark operators with nucleon spin. Intro- 
ducing dimension five operators for left-handed doublet 
ipQ and right-handed singlets ip u and ipd, members of the 
first generation of quarks 



C 



abc 



t—Q,u : d 

at the nucleon level we have 



^2 Vilpijadbdclpi, 



(12) 



771, at = a u (i] u + 7]q) + a d (r)d + T]q) (13) 
V2.N = b u (Vu - Vq) + b d(.Vd - Vq) + b -t£, 

where ??i(2), n are the rji and 772 couplings for nucleons 
defined in (7). Note that £ enters only in the 772 coupling 
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for nucleons because both are even under the charge con- 
jugation. In (13) a Ut d and b u _d are the matrix elements 
that could be obtained as the moments of the experimen- 
tally measured structure functions: a d ~ 0.4, a u <~ 0.1, 
bd ~ 0.1, b u ~ —0.05 for the neutron and charge inverted 
values for the proton. To relate the photon operator with 
nucleon, we use the simplest vector dominance model and 
obtain at one-loop level 6 7 ~ 0.13a/(47r) for neutron and 
6 7 ~ 0.24a/(47r) for proton. This results in the limit: 

\{m - vq) - oMvu - vq) + io 3 ei < kd 8 (w) 

Barring accidental cancellations, one can place separate 
limits on r\ u ^ d — ijq at 10~ 8 and on £ at 10~ 5 level. The 
orthogonal combinations r/ u ^ + r\Q are less constrained 
because they enter only in the quadrupole coupling be- 
tween the nuclear spin and external direction, and thus 
are suppressed by an additional factor of \rii \ <~ 10~ 3 [7]. 

So far we have neglected the fact that the low energy 
values for the couplings rji and £ taken at the normaliza- 
tion scale of 1 GeV do not coincide with the high-energy 
values for the same couplings generated at Mp\. With 
a simple one-loop analysis of the renormalization group 
equations, we find that several bounds can be strength- 
ened. Leaving the details for elsewhere, our results are: 
|»?Q,u,d|,|£l £ 10~ 6 and \r} e L R \ < 10" 5 with the same 
normalization on 1/Mpy. The constraints may also be 
improved by assuming a GUT scenario. 

In summary, we have shown that effective field the- 
ory provides a framework where one can derive stringent 
bounds on Planck scale interactions from terrestrial ex- 
periments. The resulting limits enhance and generalize 
the bounds obtained previously in [7] . They are generally 
far more sensitive [7] than those previously derived by 
considering astrophysical phenomena [3] where the typi- 
cal sensitivity is O(l) in units of Mp^. Actually the bire- 
fringence induced by £ 7 can be used to bound |£| < 10~ 4 
with astronomical observations [5], but our bound im- 
proves on this result by roughly an order of magnitude. 
Recently [4] infers a striking bound Vl^Vr > — 10~ 9 
(assuming rr£ — rj R — 0) from the observation of syn- 
chrotron radiation from astrophysical sources. Here our 
results provide an upper bound for positive couplings. 

Following [7], we compare our results to semi-classical 
calculations which have appeared using loop quantum 
gravity [11-13]. Examining the interactions induced for 
a Dirac fermion [12] shows that 771 = and that while 772 
is nonvanishing, it is suppressed an additional small pa- 
rameter related to the coherence scale of the gravitational 
wave function. Hence these calculations seem to be in ac- 
cord with the stringent bounds derived for the fermion 
operators. In contrast, the results of [11] suggest that £ 
should be 0(1), which stands in stark contradiction with 
the bound derived here. Hence our bounds seem to be 
in conflict with the present calculations [11-13]. How- 
ever, the latter must be viewed as somewhat heuristic. 
So the immediate challenge for loop quantum gravity is 
to provide rigorous predictions that may be compared 



to experiment. Despite the many unresolved theoretical 
issues, it is truly remarkable that precision experiments 
can already confront quantum gravity calculations with 
concrete observational bounds. 
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